Controlling spatiotemporal chaos via small external forces 
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The spatiotemporal chaos in the system described by a one-dimensional nonlinear drift-wave equation is controlled 
by directly adding a periodic force with appropriately chosen frequencies. By dividing the solution of the system 
into a carrier steady wave and its perturbation, we find that the controlling mechanism can be explained by a 
slaving principle. The critical controlling time for a perturbation mode increases exponentially with its wave 
number. 
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a^ ■ 

G^ ■ I. INTRODUCTION 

&X)| Spatiotemporal chaos(STC) can appear in a large variety of systems such as hydrodynamic systems, plasma devices, 
^ laser systems, chemical reactions, Josephson junction arrays and biological systems. Because of the potential applica- 
[ tions, controlling STC in these systems has been given much attention by scientists and technologists in recently years 
0- Generally speaking, there are two kinds of methods to control STC: the feedback control and the 

nonfeedback one [|2j. For the former one, it needs a fast responding feedback system that produces an external 
~^ ' signal in response to the system's dynamics. On the contrary, for the later one, the form and the amplitude of the 
^ I external control signal can be adjusted numerically or experimentally by directly observing the response of the system 
to the applied signal. The advantages and disadvantages of the two kinds of methods are: for the feedback one, the 
' reference state is always an unstable trajectory locating in the chaotic attractors, the controlling input is very small 
00 if the system can be well controlled, in this way, however, one must know some prior knowledge of the system before 
controlling; In contrast, for the nonfeedback one, it does not need any prior knowledge of the system, it is very easy 
2^ in practice, and particularly convenient for experimentalist, but the target states are not the unstable periodic orbits 
of the chaotic attractors, hence, the controlling input does not vanish when the system is under control. 



OO 



0^ 



Ch In this letter, we will give an example of controlling STC by adding a small external force on a system and discuss 



the mechanism. In fact, the methods of adding a pre-determined driving force or modifying an accessible parameter in 
a pre-determined way on the system to control chaos might be firstly used by Alekseev and Loskutov |]l3| , who studied 
O how to control a chaotic model of a water ecosystem by small periodic perturbations. After Alekseev and Loskutov, 
J3 Lima and Pettini jlj] successfully used the parametric perturbation method to eliminate chaos in DufFing-Holmes 
oscillator. Some experimental works about the method have been done by Azevedo jl^,Fronzoni Ding jl^] 



and Ciofini |18|. Moreover, Braiman and Goldhirsch |ig| directly added a weak external periodic force to Josephson 
. 5^ I junction system to suppress chaos even in nonresonant circumstance. Recently, the method has been intensively and 
• extensively used to control chaos not only for nonautonomous systems but also for autonomous systems and iterated 
?H maps However, within our knowledge, the method has not ever been used to control STC. This letter will give 
- - 1 the first example to use the method to control STC in a partial differential equation systems, and find for the first 
time that the mechanism can be explained by the slaving principle. 

The letter is organized as follows. First, we will briefly introduce the driven-damped nonlinear drift-wave equation 
which we used as the model equation. Then, the control method and the numerical simulation results will be given. 
Next, we will analyze the controlling mechanism. Finally, a discussion and conclusion will be given. 



II. MODEL 



The model we used in this letter is a one-dimensional nonlinear drift-wave equation driven by a sinusoidal wave. It 
reads [|| 

d<t) d^(t> d(f> d(j) M ■ t n\ 

where the 27r-periodic boundary condition, 0(a; -f- 27r,t) = (j){x,t)^ is applied. In this study, we fix the parameters 
a = —0.2871, 7 = 0.1, c = 1.0 and / = —6.0 , only the forcing strength e and phase speed Vl of the sinucoidal wave 
are the control parameters. Without the dissipation and driving terms(7 = 0, e = 0), Eq. (|l|) describes the nonlinear 
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drift-wave in magnetized plasmas After introducing the driving and dissipation terms, the competition among 
dispersion, dissipation, driving and nonlinearity makes the system display rich dynamic phenomena pl| ] p^ . 
In plasma physics, an integral quantity E{t), which is well known as the "energy" of the system, 

Eit) = l-l'^lm.,t)--aif^)^]d^, (2) 

is conveniently used to monitor the dynamics of the system. E{t) is a positive constant at e = 7 = and monotonously 
decreases to zero by the dissipation if e = and 7 > 0. 

In the present paper, we focus on STC states. On the parameter plane e — f2, the dynamic behavior of the system 
has been extensively studied by one of us (He) ||2l] , ^ . According to Fig. 6(b) in Ref. [g3|, one can find that in certain 
regions of fl, when the driving e is strong enough the system displays a STC state. It is necessary to point out 
that so called space-time chaos is an ambiguous terminology. There is evidence to show that the spatial incoherence 
may be caused by different physical effects. In previous work p^ , Hu and He reported an example of controlling 
one type of STC states by feedback method, its erratic spatial behavior is due to the overlapping of the regimes of 
the Hopf instability of different modes. In the present work, we will focus our attention on a STC state, e.g. at 
e = 0.22, ^2 = 0.65, the transition to STC is caused by a crisis |p3|| . 



III. METHODS AND NUMERICAL SIMULATION RESULTS 



For controlling STC in the system, we directly add a small temporally periodic signal on the right hand side of Eq. 
(I). Thus Eq. (§) is changed to 

d4> ^ 9cp del, w o.^ ^ ^ ,^ 

m + + + ^^9^ = -^'^ - + VCOs{u;t), (3) 

where rj and to are the control parameters. 

The pseudospectral method |Q is used to simulate Eq. (||). In the numerical simulation, we divide the space into 
128 points. 256 and 512 points have also been used, the results do not show qualitative difference. Therefore, the 128 
points are used in our work. 

In order to find suitable rj and uj to control STC, we first calculate the power spectrum of (j){xo,t), denoted by S{uj), 
at a fixed point x = Xq in space. From the power spectrum S{uj) one can find the characteristic frequency region of 
the chaotic attractor. For example. Fig. 1(a) shows the spectrum of the STC state to be controlled, from this figure 
one can see that except for the frequency of driving wave(17 = 0.65), there are many peaks in the region oj £ [0.2, 0.9]. 
Our experience shows that the appropriate u for a successful control should be chosen in this region. Now, fix an w 
in this region, one can try to find a suitable rj. We find some windows in w — 77 plane, with which parameters the 
turbulent motion of the STC state can be successfully suppressed. The results are given in Fig. 2, in which the circles 
stand for the suitable parameter points for controlling STC in system (1), while the crosses indicate the parameter 
points at which the turbulence are failed to be controlled. In the following, we choose uj — 0.756 and rj = 0.1 as an 
example. To compare the dynamic behavior of the system before and after the control, the power spectrum, energy 
and spatiotemporal patterns are calculated. Figure 1(b) shows the power spectrum S{ui) at a fixed space point (a; = tt) 
after the control. In contrast to the wide noisy spectrum in Fig. 1(a), here only a few lines are left. In Fig. 3, We 
give the time series of the energy E(t), one can see it is chaotic before control, which becomes harmonic oscillations 
after a while when the control signal is added. The spatiotemporal patterns of the system before and after control 
are shown in Fig. 4(a) and Fig. 4(b), respectively, obviously the turbulence is strongly suppressed . 



IV. MECHANISM AND SCALING BEHAVIOR 

To analyze the mechanism of the control method introduced above, firstly, let us observe the solution of Eq. (^ in 
the frame of the driving wave. By setting z — x ~ fit and t = t , Eq. (^) can be rewritten as 

9 52 Q Q2 QM Q± 

+ ^-Sr?)^ = + °'^^)'t' ^ esinz rycos(wr). (4) 

OT oz'^ oz oz^ oz oz 

Then, we divide a solution of Eq. ^ into a carrier steady wave 4>o{z), and it's perturbation S(t){z,t), i.e.. 
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^{z,t) =(f>o{^) + ^Hz,t). (5) 
Here (po satisfies the steady equation 90o/9r — 0, i.e., 

o ^ n , \j, ^'^^o r, 9(1)0 . 

52— (1 + a-—^)(po - c— /0o^ 700 - esmz = 0. (6) 

oz oz^ oz oz 

Insert Eq. (||) into Eq. (^), and using Eq. one has 



1 + ""H^ '^'^ = 1 + "^^'^^ ~ f^o-^ /O'/'^ fScf) „ 

OT oz"^ oz az^ oz oz oz oz 



—jS(f> + ricos{u!T). (7) 
By making the Fourier transformation of (I)q{z) and S(j){z,t) as 



(j)oiz) ^ Ao + liniN^oo Y.k=i cos(fcz + 9k), 

6(j){z, t) = ba{t) + liniN^oo Ylk=i ^k{t) cos(fcz + Qfc(i)), 



(8) 



Eq. (||) is transformed to a set of infinite algebraic equations, and Eq. (|7|) a set of infinite ordinary differential 
equations 

= -T^'o + JycosM), 

^ = Nkit) + fkboAk siniOk - afc)/(l - afc^), (g) 

dofc _ llf (-j-N _ /fc(Ao+bo) _ fkbpAk cos(0k-ak) 
— ^'>'^k\i') (l-afc2) (l-afe2)6fc ' 

where 

, fk 



Nk{t) 



l-aP 2(l-afc2) 
X { ^ [Aihv sin(6li + a;/ - au) + ^kbi, sm{ai + av - ak)] 

l+l'=k 

+ ^ [Aibi> sm{9i - ai> - ak) + ^kbu sm{ai - av - ak)] 

l-l' = k 

+ ^ [Aibi^ sm{-9i+ av ~ ak) + ^bibii s\-a.{-ai + av - ak)]}, 



l'-l = k 



Mkit)^ ~ki-^-n) 



l-aP ' 2{l -ak^)bk 
{ ^ [Aibv cos{9i + av - ak) + ^kbv cos{ai + av - ak)] 



X 

l+l'=k 



+ ^ [Aibv cos{9i - av - ak) + ^bik' cos{ai - av - ak)] 

l-l'=k 

+ ^ [Aibv cos{~9i + av - ak) + ^kbv coa{-ai + av - ak)]}, (fc = 1, 2, . . . , TV). 

l'-l=k 

In numerical simulations, the mode equations for Ak,9k and bk,ak can be solved simultaneously with truncation. 
The proper N should satisfy the condition that the modes k > N give negligible contributions to the actual solution of 
the equation For the present parameters we use N = 13, i.e., the 27-dimcnsional ordinary differential equations 
will be solved. With these mode equations one can obtain qualitatively the same results as the direct simulating Eq. 

Now let us analyze the numerical results of the solution of Eqs. (W. From the time evolution curves of bk, one 
can see that, for every mode k, there exists a turning point r^, before which b^ shows chaotic motion with larger 
amplitude, while after which it becomes a quasi-periodic vibration with much smaller amplitude. As an example. 
Fig. 5(a) depicts the time evolution of 63 and bn in the beginning 50 time units after the control signal is added 
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and another 50 time units when they behave quasi-periodically. In order to determine , let us calculate the relative 
average amplitude value of each mode 

rk{t)^mt)~bk)/bk, (10) 

where bk is the time average value of bk in the quasi-periodic states, &^ is the points only on the Poincare section 
pGf . As an example, the results of rs{t),rio(t) and ri2{t) are shown in Fig. 5(b), one can see that, after adding the 
periodic control signal, rk(t) starts to decrease exponentially, but they never tend to zero. After a certain time point, 
the oscillation level of each rt no longer decrease. 

Next, we will study the dependence of the critical point on the mode number k. In the decreasing phase of r^, 
its maxima decreasing linearly in the logarithmic diagram, while in the quasiperiodic phase, the maxima of bk are 
about a constant. The intersection point of the two straight lines gives t^. In Fig. 5(c), we plot as a function of 
mode number k in the semilogarithmic plane, the points can be well fitted by a straight line. Consequently we obtain 
a scaling law 

r*k oc eP\ (11) 

here by the least square method we get (3 = 0.1169. When amplitude bkS transit to the quasi-periodic phase, the 
mode phase ,afe, of the originally chaotic mode also become very quiescent. We can also determine the turning points 
from the evolution of ak^ the same scaling law as Eq. (pi] ) is obtained, with j3 = 0.1158. 

From these results we conclude that, by using our time-periodic signal to control a STC state, the modes{6fc, a^} 
will be controlled one by one from the smallest wave number k to the highest one. The critical time for a mode to 
be controlled increases exponentially with the wave number k. This phenomenon can be understood by the slaving 
principle in synergetics [0. From the first equation of Eqs. (^, one can see that, in fact, the external force is added 
on the mode fc = 0, the evolution of 69 is uniquely determined by the applied control signal and the dissipation rate 
7. Its motion does not influenced by the other k ^ Q modes. The motion of 60 is adjusted to periodic one very 
quickly after the control signal is added. Then the other modes fc 7^ are slaved by b^. This fact can be seen from 
Eqs.(^ for fc 7^ 0. In both equations for bk and aj,, there is a term depending on b^, which is the contribution of the 
mode-mode couplings between b^^Ak.Ok} and {bk;ak}- Through these terms 60 plays the slaving role. Meanwhile, 
since the modes of the steady wave 4>q{z), {Ak, 9k} are constants in the frame of {z, r}. the modes {bk, ak} also play 
a significant role in the solving process. Moreover, the numerical results tell us that bo is a slowly relaxing variable 
and the other modes are fast ones in the beginning stage when the control signal is just added. Thus, bo plays an 
important role as an order parameter in the system, and it is affected by the control force. When the control signal is 
added, 61 and ai are immediately slaved by bo- Since the effect of bo on {6^, ak} is through {Ak, Ok} and Ak decrease 
with k ^ 00 , the slaving intensity is not unform for different modes. From Eqs. (||), one can see that the smaller the 
mode k, the higher the slaving intensity it has. So the mode fc = 1 can be controlled firstly. Then, it starts to slave 
the modes A: = 2,3, until all the modes are tamed. In the whole process, the system changes from an extremely 
disordered(a STC state) into an ordered one by the self-organization principle. 



V. CONCLUSION AND DISCUSSION 



In conclusion, we have successfully suppressed very turbulent STC states by directly adding a temporally periodic 
force on the system. The appropriate input frequency can be chosen by analyzing the spectrum of the turbulence. 
In contrast to the feedback method for which one needs the knowledge of the target state, the present method is 
easy to be realized in practical experiments. The controlling mechanism is studied by investigating the evolution of 
perturbation wave to the carrier steady wave. We find that the perturbation modes are controlled one-by-one from 
long to short wave length by a slaving principle, and the critical controlling time for the mode increases exponentially 
with its wave number k. Since the slaving effect is very important in nonlinear science, these results may be interested 
in by scientists and engineers. The schemes of how to select the controlling frequency may also be helpful for the 
experimentalists . 
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FIG. 1. Power spectrum of potential at fixed z{z — tt) before(a) and after(b) control. The time interval,Af, between 
consecutive samples is 0.5, The 1024 points are used to do the FFT, the total sample points are 4096. 
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FIG. 2. Global view of control parameters on the plane cu — rj. Points o and x stand for the suitable and unsuitable 
parameter points for controlling STC in system (1), respectively. 
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FIG. 3. The energy E(t) as a function of time before and after control. 
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FIG. 4. Spatiotemporal patterns of 4>{z,t) before(a) and after(b) control. 
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FIG. 5. (a) The evolution of the amphtudes, bk{t){k = 3, 11), for the perturbation wave, the dashed hnes scaled by the 
right axises indicate the external force, (b) The evolution of rk{k = 8, 10, 12). (c) The scaling law of r^. 
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